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ABSTRACT 


This paper presents electromagnetic data recorded from lightning strikes. 
The data analysis reveals general characteristics of fast electromagnetic fields 
measured at the ground including rise times, amplitudes, and time patterns. A 
look at the electromagnetic structure of lightning shows that the shortest rise 
times in the vicinity of 30 ns are associated with leader streamers. Lightning 
location is based on electromagnetic field characteristics and is compared to 
a nearby sky camera. The fields from both leaders and return strokes have been 
measured and are discussed. 

The data were obtained during 1978 and 1979 from lightning strikes occur- 
ring within 5 kilometers of an underground metal instrumentation room located 
on South Baldy peak near Langmuir Laboratory, New Mexico. The computer con- 
trolled instrumentation consisted of sensors previously used for measuring the 
nuclear electromagnetic pulse (EMP) and analog-digital recorders with 10 ns 
sampling, 256 levels of resolution, and 2 kilobytes of internal memory. 


I . INTRODUCTION 


In the summers of 1978 and 1979 a set of measurements of the electromag- 
netic fields from natural lightning was made near the top of South Baldy peak 
(elevation 3275 m) which also houses Langmuir Laboratory near Socorro, New 
Mexico, U.S.A. As will be discussed in this paper the data consist of time 
waveforms of about 20 ys duration at a 10 ns sample spacing for the three elec- 
tromagnetic field components (one electric and two magnetic) which exist 
adjacent to a conducting ground plane on top of the "iron Kiva M (or Kiva) (our 
screen room buried in the top of the mountain) and the surrounding earth surface 
(see fig. 1.1). The details of the instrumentation system are discussed in 
another paper (ref. 5). 

As background, these measurements resulted from some interest in the dif- 
ference between lightning and the nuclear electromagnetic pulse. For nuclear 
detection purposes it is desirable to reliably detect a nuclear EMP event 
without false triggering by lightning. The ARGUS-1A detector system was 
designed with the high altitude EMP environment in mind (ref. 2). For this 
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purpose magnetic (loop) sensors were used to avoid the low-frequency electric 
field from lightning and local static charging effects. In addition the sen- 
sors respond to the derivative of the field to emphasize the high frequencies 
present in the high-altitude nuclear EMP. Essentially one is looking for a 
pulse with 10 ns to 100 ns important characteristic times. 

In order to understand the sensitivity of the ARGUS-1A to lightning two 
things were done. First the detectors were placed on South Baldy peak and their 
response to lightning monitored. Second, and more important for this paper, the 
fast transient characteristics of the transient electromagnetic fields (in time- 
derivative form) were measured. It is these measurements and some analysis of 
them which are presented here. 

The field sensors consist of two MGL-3 B-dot sensors and one ACD-5 D-dot 
sensor located on the roof (ground plane) of a buried shield enclosure (Kiva) . 
The roof of the Kiva is flush with the earth surface and is electrically con- 
nected to the center of a 30 m x 30 m wire-mesh ground plane. Data recording 
is accomplished by Biomation 8100 waveform recorders controlled by an HP 9825 
calculator. A more complete description of the instrumentation is presented in 
reference 5. 

At first we had little idea of what to expect. We made some measurements 
(section 2) which showed us some of the signal levels and patterns in the wave- 
forms of the field time derivatives, and in the time integral of these wave- 
forms. Based on the kinds of pulses observed and their time patterns it 
appeared that these related to individual fast pulses along some leader forma- 
tion process. This led to a model (section 3) which with ranging information 
(section 4) allowed us to infer some of the characteristics of the current in 
the lightning streamer. Section 5 then applies this to some of the data 
obtained . 


II. GENERAL CHARACTERISTICS OF OBSERVED ELECTROMAGNETIC FIELDS 


With the instrumentation at the Kiva discussed in reference 5 measurements 
were begun in the summer of 1978. For use in our data analysis the 3$/8t 
sensors (MGL-3) had A^ eC j = 0.1 m^ , and the 3l5/3t sensor (ACD-5) had A ee q= 1.0 m 
for use in the formulas 


v h = \ 


eq 


at B 


(open circuit voltage) 


I = A • — D (short circuit current) 
e e 3t 

eq 


( 2 . 1 ) 


where the fields are those in the presence of the ground plane (i.e., including 
reflection). The response times (10-90 rise times of integrated output for 
step-function exciting field) into 50 loads are short compared to our 10 ns 
sampling resolution and hence neglected. 
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The three sensors responded to the north and east components of 3$/3t and 
to the vertical (up) component of 3D/3t. The signal from 3B/3t north was atten- 
uated a factor of 2 by a power splitter used to send this signal to an addi- 
tional recorder for other purposes. These corrections and the equivalent^ areas 
are all removed for the displays of 3B/3t and 3^/3t in this paper. The B and D 
waveforms are determined by numerical integration of the digital data provided 
by the Biomation 8100 waveform recorders. Note that the integrated waveforms 
exhibit baseline slope which is at least partly due to errors in determining 
the baseline values for the original waveforms; this baseline shift is partly 
removed by estimating the baseline position near the beginning of the recording. 
The long-time variation of the integrated waveforms should not be trusted, 
nevertheless; the purpose of the integration is to see the shape and amplitude 
of the fast pulses which are more reliable. 

A nice property of the Biomation 8100 waveform recorder is that one has a 
significant length of recording time before the pulse which triggers the system 
is recorded. This pretrigger data time was set at about 4 ys. The self 
trigger was used on all three channels with the first trigger level being 
exceeded triggering all the recorders. Triggering in this case means freezing 
the data stream and transferring the 20 ys of information to magnetic disk. 

Since the recording is digital there is also some limitation in the signal 
amplitude (voltage) resolution. The full scale is 256 uniformly spaced levels 
with zero approximately centered. Typical signals deviated about one third 
from baseline (zero) to the full level on one side of zero. 

For some preliminary information let us now consider two examples of the 
kinds of waveforms we have seen. Detailed analysis of such waveforms is con- 
sidered in section 5. A representative example of a leader waveform is given 
in figures 2.1 and 2.2. A representative example of a return-stroke waveform 
is given in figures 2.3 and 2.4. In both cases only the east components of 
3$/3t and of S are exhibited. 

Consider the leader waveform in figure 2.1. Note the complex structure 
with many individual pulses indicating the progression of the leader as a 
sequence of impulsive events. The horizontal distance (by flash to bang) to 
this stroke was about 350 m. 


Note the large pulse at about 4 ys which triggered the recording. The 
first 7 ys of the record is expanded in figure 2.2. From a printout of the 
digital data this pulse has the approximate characteristics: 


pulse at 4.0 ys 

peak - 2.42 T/s 
zero to peak rise ^ 90 ns 
10-90 rise - 50 ns 

( 2 . 2 ) 
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B 


E 


! 


baseline to peak - 0.169 yT 


baseline to peak rise - 150 ns 


width - 500 ns 


10-90 rise - 80 ns 


T - 70 ns 

Here T has been defined as a characteristic time for the rise portion of the 
waveform with the formula 

__ waveform peak ^2 

derivative waveform peak 

Selecting a few more pulses from this waveform for comparison we have 
pulse at 1.3 ys 



(2.4) 


B. 


E 


! 


baseline to peak ^ .079 yT 


baseline to peak rise ^ 100 ns 


width ^ 250 ns 


10-90 rise ^ 60 ns 


X - 53 ns 


and another 


pulse at 12.9 ys 



(2.5) 
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{ baseline to peak - -.03 yT 
baseline to peak rise - 80 ns 
10-90 rise - 40 ns 
width - 160 ns 

T - 35 ns 


and another 

pulse at 



14.0 ys 

peak - -.547 T/s 

zero to peak rise - 40 ns 

10-90 rise - 30 ns 


! baseline to peak - -.014 yT 
baseline to peak rise - 50 ns 
10-90 rise - 35 ns 
width - 70 ns 


( 2 . 6 ) 


T - 26 ns 

Comparing these individual pulses we have seen characteristic times for the 
rise ranging from about 70 ns to less than 30 ns. Also one can note that there 
seems to be a general trend in that smaller pulses have shorter characteristic 
times for the rise and shorter pulse widths. Pulse widths vary from about 
500 ns to about 70 ns. 

Now consider the return-stroke waveform in figure 2.3. The 3B/8t structure 
is somewhat simplified by comparison to the previous leader example. The hori- 
zontal distance (by flash to bang) to this stroke was about 850 m. Comparing 
the peak derivative signals between return stroke (fig. 2.3) and leader (fig. 
2.1) note that they are of comparable magnitude. Assuming that the leader is 
elevated somewhat above the ground, the inferred slant ranges to the sources 
for the two examples are crudely of the same order as well. It therefore seems 
that we may typically expect peak derivative waveforms from the larger leader 
pulses to have about the same strength as those for return strokes (at compar- 
able distances) . 

The main pulse at about 4 ys which triggered the recording is expanded in 
figure 2.4 for the first 7 ys. Note the wide pulse characteristic of the 
return stroke as indicated in the east component of the H waveform. From a 
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printout of the digital data this pulse has the approximate characteristics: 


pulse at 4.0 ps 

{ peak - -2.81 T/s 

zero to peak rise - 110 ns 
10-90 rise - 35 ns 

(2.7) 

! baseline to peak - -.129 yT 
baseline to peak rise - 160 ns 
10-90 rise - 60 ns 
width - 1000 ns 


T - 46 ns 

It does seem that return strokes can be somewhat fast. We have looked at other 
return strokes and noticed similarly small values of T. For this type of pulse 
of B the concept of width is somewhat problematical because of the current that 
continues in the arc after the fast part of the pulse has occurred. 


III. FAST TRANSIENT ELECTROMAGNETIC FIELDS AT EARTH SURFACE 
RELATED TO INDIVIDUAL LIGHTNING CURRENT PULSES 


Having seen some of the transient waveforms in section 2, and noting that 
many of the individual pulses are extremely fast with rise times of the order 
of some tens of ns, it appeared that one might be able to model these pulses as 
events somewhat localized in space. This concept uses the limitation of the 
speed of light on current propagation as compared to the distances of hundreds 
to thousands of meters from source to observer. As will appear later, this 
approach yields some valuable insight into the lightning streamer currents and 
propagation. 

Let us begin by defining a spherical (r,0,cf)) coordinate system based on 
the Kiva as in figure 3.1A. For our purposes (r,0,(f>) signifies some position 
in space which is the source of the observed lightning signals at (near) 
r = 0 (the Kiva) . In addition to the spherical coordinates we have cylindrical 
(T,(j),z) and cartesian (x,y,z) coordinates all related via 
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x = T cos(cj)) = r sin(0) cos((j)) 
y = ¥ sin(cf)) = r sin(0) sin(cj)) 
z = r cos(0) 


(3.1) 


As indicated 


2 u / 2 x 2 

r = i + z 

2 2 2 

¥ « x + y 


sin((j)) _ cos(c})) 
Y x 


with x and y taking the signs of cos(c()) and sin(cf)) respectively, 
in figure 3.1A the coordinate system is oriented so that 


x is north (geographical) 
y is west 
z is up 

0 is zenith angle 

(J) is counterclockwise from north 


(3.2) 


Associated with these coordinate sets there are sets of unit vectors 
(right-handed) as 


1X1=1, 1x1=1, lxl=l 
xyzyzxzxy 


4x1=1 , 1 . x 1 =11 X 1 = 1 

r <p z <p z r z r <f> 


1 X 1 = 1 1 X 1 = 1 , 1 X 1 = l , 

r 0 <J> * 0 cj) r (f) r i 


(3.3) 


which are related by 


ly = l x CO8(0) + 1 sin (cj)) 

1, = -1 sin(cf)) + 1 cos(d)) 
<p x y 


l=l m sin(0) + 1 cos ( 0 ) 
r 1 z 


(3.4) 


lg = cos(0) - 1^ sin(0) 


In addition to the coordinates for the source point with respect to the 
observation point let us define a set of unit vectors for describing the wave 
propagation to the observer. These are a right-handed system 


1 1 X 1 2 = X 3 ’ 1 2 X 1 3 = h * l 3 X l l = 1 2 


(3.5) 


which are related to our previously defined coordinate unit vectors as 
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( 3 . 6 ) 



-* -* 



Viewed from an observer at the Kiva one has the picture in figure 3 . IB in which 
->• 

1^ = direction of propagation 
= pointing to observer 
->■ 

I2 = horizontal to right ( 3 . 7 ) 

1^ = "up" as modified by being perpendicular to 1^ 

*► — 

Note that I2 and I3 are polarization vectors for the wave. The plane of inci- 
dence is defined by^r and the z axis; it contains Tj and I3 as well as l r , *1 q , 
Ty, and "T z ; l 2 and I9 are perpendicular to the plane of incidence. 

Next consider some distribution of current density J(r,t) in free space. 
The electromagnetic fields from this can be represented as (ref. 1 ) 


E. 

inc 


H. 


inc 


(r T , s) = - <s^Z (r , r" ; s) \ J(r", 
( r \ 8 ) = ~<V xt(r\r"; s) ; 





( 3 . 8 ) 


with the subscript inc indicating incident at the observer at r ! = 0 and where 
r^r" are now general coordinates (positions), s is the Laplace-transf orm vari- 
able (or complex frequency) with respect to time t (in general two sided) and 
the tilde - indicates a Laplace-transformed quantity. The impedance kernel is 


Z(r ,r M ;s) = sy G (r 1 ,r";s) 
o o 

<f (r',r";s) = T - Y^V'V'l G (r\r";s) 

O L Jo 

= dyadic Green T s function of free space 

. ( 3 . 9 ) 

G = Y 7~ (for r f ^ r") 

E scalar Green’s function of free space 

g 

y = — = free space propagation constant 
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c = 


3 x 10 8 tn/s 


V\x z 

o o 


= speed of light (in free space) 


Z E /— - 377 Q 
° J £ o 


= wave impedance of free space 

. ->■ •> . 

C = y I r T - r" | 

The symbol <,> indicates integration over the space coordinates common to the 
two terms (r ,f in (3.8)). The free-space dyadic Green’s function has the 
explicit form 

t(K> ; s) . x e- C j[-2 C - 2 - 2?' 3 ]y R 

+ [ c - 3 + c - 2 + 5 - 1 ][T-T R T R ]j + -i I 6( 


r - r M ) 1 


(3.10) 


1 = identity dyad 

-r T — r- n 

I E— — forrWr" 

r T - r" I 


where the delta function is introduced so that the remaining part can be 
handled in (3.8) as a principal value integral with a small spherical zone of 
exclusion centered on r T = r" (ref. 3). 

For our present purposes we assume that the distance from source to 
observer is large compared to the size of the source, i.e., 


. -> 
r = 0 


r" = |r ,f | 


(3.11) 


where r is taken as some effective center of the source. Furthermore we assume 
that frequencies of interest are sufficiently high that only the leading term 
for large r ,f (the |r f - r u |“l or term in (3.10) need be included giving 


^ + e ~Y I r ? " r"| 

G (r',r”;s) - — 

o . i-> , 


4tt I r T - r"| 

Y \- r " 

— — i 

4irr t 


(3.12) 
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(3.13) 


1 = 1 - l n l_. = transverse dyad 

t K R 


Now referring to figure 3.1 we can approximate as 
-> 

1*1,= "I 
R 1 r 

-«-► -y «->■ -> -* 

1 * 1 - I! = 1 - l 1 

t 11 r r 

This allows us to write (from (3.8)) 

-Yr 

- ye' 

E. (0,s) = - 1 • T(s) 

inc 4?rr t 


?(s) = s 


Y(r-1 *r") ~ ^ 
e r J(r",s) dV" 


In time domain we have 

3 


?(t) 


3t 


V 


J(r",t + -(r - t • r")) dV" 
c r 


(3.14) 


(3.15) 


This T can be thought of as a kind of effective source vector. Assuming that 
changes in J propagate over the source volume at velocities slow compared to c 
we have 


T(s) 



J(r" , s) dV' 



J(r M 


V 


t) dV n 


(3.16) 


The magnetic field incident at r T = 6 can be found by expanding from (3.8), or 
by simply noting that the far field (|r T - r"|~l term) has a plane-wave rela- 
tion between the electric and magnetic field as 


H. (0 , s) - ~ t, x E. (0,s) 
me Z 1 me 

o 


-yr 

e 1 -y -y 

7 1, x T(s) 

47rrc 1 


In time domain the fields are 


-y 

E. 

me 


(0,t) * - 


4Trr t 



(3.17) 


(3.18) 
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H. (0,t) - - yi— t. x f(t - -) 
inc 4irrc 1 c 


Near r T = 0 the incident field is an approximate plane wave propagating in 
the "t]_ direction; it can be represented as 


1 . r ' 1 • r' 

E. (r ' , t) - E. (t ) 1 + E (t ) 1. 

me h c 2 e c 3 


H. (r\t) = -t Je, (t 

me Z / h 

o ( 


1 • r 9 1 • 

) t - E (t 

c ' 3 e c 


-) t 


(3.19) 


This decomposes the fields incident at the observer into two parts, an E wave 
(or TM wave) denoted by a subscript e, and an H wave (or TE wave) denoted by a 
subscript h. This decomposition is based on the polarization with respect to 
the z axis (i.e., a TM wave has the magnetic field perpendicular to the z axis, 
and a TE wave has the electric field perpendicular to the z axis) . 

Now when the incident fields reach the ground plane (the x,y plane), the 
conducting ground plane reflects the fields so that (at least at the higher 
frequencies because of the finite dimension of the ground plane on the earth) 
approximately we have 


E - 

21 

-> 

• E. 



z 

z 

me 




->• 

->■ 



H - 

21 

• H. 

= K = 

-H 

X 

X 

me 

N 



-> 

-+ 



H - 

21 

• H. 

= IT, = 

-H. 

y 

y 

me 

w 



(3.20) 


The remaining three fields components are approximately zero. Note the orienta- 
tion of the cartesian coordinates with respect to the (geographical) compass 
positions (north = N, east E E, south E S, west E W) . 

On the ground plane we then have 

E = 2 sin(0) E 
z e 


% - - H 2 = T E e 

o 


(3.21) 


Hy = -cos(6) H 3 = - — cos(0) E h 

o 


Then construct 


tE = sin(0) E 
/ z e 


49 


(3.22) 


Z Z 

- 7 T- H , = {-H sin(<J>) + H cos((}))} = E 

2 <p 2 x y e 

Z Z 

^ Hy = {H x cos(c J>) + sin(<()) } = -cos( 0 ) 

From our measurements we have three pieces of information: E , H , and H . 

We have four unknowns: E , E, , 0, and <b . x Y 

e h 

If we knew the source angles ( 0 , cf) ) , or at least one of these, then the two 
waves E e and E^ could be constructed. Let us for the moment assume that 0 and 
(f> are known. From (3.18) and (3.19) together with E e and E^ we can find 



4iTr + 

y inc 
o 


( 0 ,t) 


47Tr | 
1 


E h (t) T, 


+ E e (t) T 3 | 


(3.23) 


This shows that we can obtain the relative strength of the transverse (2 and 3) 
components of the source T by comparing E^ and E e . Referring to figure 3. IB 
one can plot on such a diagram an angle that 'f makes with respect to say ^3 as 
another vector on the plane perpendicular to T]^. Note that the 1 component of 
is not obtained in this technique since it does not appear in the far-field 
expansion at r 1 = If one knew r then the 2 and 3 components of t 5 would be 

completely known. Techniques for estimating r are discussed in another section 
of this paper. Having r then one can use (3.23) to find T 2 and T 3 ; this is 
used for later plots of the transverse part of T. 

- $■ > -> 

Having l t • T one can also make some estimates concerning currents in the 
fast pulses in the source region near r by writing 

J(r",t) dV" (3.24) 

V" ) 

Here A indicates the change in the quantity; for later plots this is the change 
from the ambient value befor^ a fast pulse to the peak of the pulse. (This 
convention is employed for 3T/3t also.) An effective streamer speed v ^ can 
be used to obtain some indication of the current I, or at least some rough 
bound on it. 

Now consider the determination of 0,<j). Figure 3.2 shows a type of display 
of coordinates which is useful for our purposes. Superimposed on the cartesian 
directions and compass coordinates we have a polar plot in which 0 is the radial 
coordinate with 0° <L 0 <_ 90°. The azimuthal coordinate is <}> with 0 <. < 360°. 

Returning to (3.22) our three measurements (E z , H x , Hy) can be used to find 
a relationship between 0 and <J> . Specifically the first two of (3.22) can be 
used to eliminate E e giving 


T 5 




= A 


IJL 

|3t 
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(3.25) 


Z 

o 


-H sin(<|>) + H 
x y 


COS((j)) 


E 

z 

sin(0) 


as our 0 , cf> relation. Constraining 0 <. 0 <. 90° and <() real gives a curve in the 
0,cf) ,f plane l! of acceptable solutions of (3.25). For each pulse we can consider 
the change in the fields (ambient to peak) in (3.25) or some other aspect of 
the pulse waveforms. Besides the field change we have used the change in the 
time derivatives of the fields (3E z /3t, 3H x /3t, 3H z /3t) in (3.25) to determine 
the 0,0 curves. Both the field changes and changes in field time derivatives 
(on the rise of the pulses) typically give comparable results, as will be seen 
later. 


A special case is encountered if the incident field is a pure H wave (i.e., 
E parallel to the ground plane). In this case E e = 0. From (3.22) we have 

H 

tan(<p) = ^ (3.26) 

X 

which has two solutions for 0 separated by 180°. In the 0,0 "plane” in figure 
3.2 this gives a locus of solutions a "straight line" passing through the 
zenith (0=0). 

Having a contour in the 0,0 "plane" resulting from a particular measured 
fast electromagnetic-field pulse (or its time derivative) one would like to 
find some point along that contour that represents the true angles 0,(j) to the 
source. One way to do this uses an assumption that another pulse in a particu- 
lar recording (about 20 ys wide in the presented data) comes from approximately 
the same source location. This is physically reasonable if one considers a 
leader propagating through the air to establish the lightning channel. If each 
pulse represents some brief extension of the arc channel, and if the distance 
the streamer propagates during the pulse is small compared to the distance to 
the observer, then the change in (0,0) is small. Furthermore, if the two pulses 
under consideration have different polarizations ^different proportions of E 
and H waves) due to different directions of l t • T, then (3.25) will give in 
general different 0,cf) contours for the two pulses. The intersection of the two 
contours can then be taken as an experimental determination of (0,(f>) for both 
pulses . 

Going further, it will be observed in considering the data that the inter- 
section of 0 ,(f) contours is remarkably successful in determining ( 0 , <f>) values. 

In some cases all the 0,<f> contours (of the order of ten contours) can all pass 
quite near a common (0,0) approximate intersection. In other cases subsets of 
the contours for a particular waveform set form a few distinct approximate (0,0) 
values. Some cases are presented in the data illustrating some of the complex- 
ity that can be unravelled in this way. 

Having determined the (0,0) values appropriate to particular pulses, and 
estimating the corresponding r, one goes back to (3.22) to find E e and E^ (or 
their derivatives with respect to time). Then (3.23) gives the transverse com- 
ponents of 'f 5 for plots with respect to the I 2 an< ^ ^3 directions as in figure 
3. IB. 
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IV. ESTIMATION OF RANGE TO LIGHTNING CURRENT PULSES 


The electric field changes produced by lightning in this study were also 
sensed with a "slow antenna" field change meter and recorded onto FM magnetic 
tape, together with the IRIG B, 1 kHz time code that was distributed to coordin- 
ate the thunderstorm and lightning measurements. This time code was synchron- 
ized with Radio Station WWV and permitted the slow measurements to be correlated 
with a time resolution of 1 ms. A Globe 100 B microphone with an extended low 
frequency response was mounted on the metal mesh adjacent to the Kiva and was 
used to detect the arrival of thunder. The thunder signals were recorded on a 
second FM magnetic tape channel. In this manner, the time of occurrence, the 
nature of the lightning flash and the time from f lash-to-thunder could be 
determined . 

Two acoustic arrays, each consisting of three microphones at the vertices 
of 30 m, equilateral triangles were located in the meadow, south of the Kiva. 
(See Figure 1.1.) The time-difference of arrival of individual thunder peals 
across each of these triangular arrays permitted a determination of the direc- 
tion cosines for the normal to the acoustic wave front. With these and the 
distance from each source to the array calculated from the speed of sound and 
the time interval from field change to peal arrival, a first estimate could be 
made for location of each source. Corrections were then made for the local 
winds and for thermal refraction in a second approximation for the source loca- 
tion using the technique described by Winn et al. (ref. 4). 

In addition to the use of these acoustic techniques for the location of 
lightning, the flashes in the vicinity of the Kiva were video-recorded with a 
time resolution of 16 ms per TV frame. One TV camera recorded the view from an 
all-sky reflecting, parabolic mirror located on top of South Baldy peak about 
50 m north of the Kiva, while another TV system located at Langmuir Laboratory 
recorded the view of lightning over South Baldy peak as seen from 2 km to the 
SSE . 


These video recordings have been played back, frame by frame, onto a video 
monitor where they were photographed together with an output displaying the 
IRIG B time code that was also recorded on the audio tracks of the video tape 
recorders . 

All of the flashes reported in this study were "cloud-to-ground" discharges. 
Our first approximation for the slant range from the Kiva to the lightning cur- 
rent element of interest is based on the assumption that the time from flash to 
thunder multiplied by the speed of sound in air gives the horizontal distance 
from the Kiva microphone to the closest part of the lightning channel (that is 
striking the earth approximately vertically) . 

From section 3, the zenith angle 0 and the azimuth angle <j> can be estimated 
for the current element from the electromagnetic signal that it produced. 

The approximate slant range r from the Kiva to the current element is 
therefore given by 
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r - (336 m/s) (At) /sin(9) 


(4.1) 


The cartesian coordinates of the current element are therefore 


x 


(336 m/s) (At) cos(cj)), toward the north 


y = (336 m/s) (At) sin(cj)), toward the west 


(4.2) 


z 


r/cos(0), above the Kiva 


In the four cases to date in which we have compared the thunder source 
maps with current element coordinates, there has been a fair agreement between 
the two somewhat independent location systems. If the range chosen from the 
horizontal distance assumption is too small, a better^estimate for the range 
can be extracted from the intersection of the vector r — directed along the 
calculated values of 9 and cf> — with the thunder channels mapped by the acoustic 
technique . 


In this section the techniques discussed in sections 4 and 5 are applied 
to some example data sets for four selected waveforms. 


Our first example is given in figures 5.1 through 5.8. It is a good 
example of leader-like behavior and is associated with a triggered lightning 
event. The lightning stroke was triggered by a wire-trailing rocket launched 
from the position indicated in figure 1,1 which is 350 m in horizontal range at 
an azimuth of <J> = 155° with respect to the Kiva. The rocket was launched 
approximately vertically. Based on visual observations its height was estimated 
as about 600 m above the Kiva when it initiated a lightning discharge. This 
corresponds to a zenith angle of 9 - 30° with respect to the Kiva. Based on 
the extreme electric-field enhancement at the rocket, it was expected that the 
streamer direction of interest was upward, beginning from the rocket. This 
polarity is assumed later in the data analysis. 

Looking at figure 5.1 we see the leader-like behavior on the 20 ys wave- 
form for the three derivative field components. For convenience 9 pulses are 
identified by numbers corresponding to their peak times of occurrence (within 
the record event) in microseconds; the north component of 9B/3t had the gener- 
ally largest pulses and was used for this labelling. Figure 5.2 results from 
the numerical integration of the digital data from which figure 5.1 was pro- 
duced. Here the pulses are also identified (on the east component of B in this 
case). Note the general drift in the baseline which may be partially due to a 
little offset in the recording of the derivative data, but probably is a result 
of a continuing current flowing up the channel initiated by the rocket. 


V. ANALYZED WAVEFORMS 


A. Rocket-Triggered Lightning 
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Figure 5.3 shows the slow changes in the electric field, first from the 
rising of the grounded wire, followed by the emission of upward-going positive 
streamers. Also shown is a recording of the narrowband RF power received cen- 
tered on 34 MHz, as measured by a different experiment (Hayenga) with a complex 
antenna pattern. For our purposes it is interesting to note that the RF power 
emissions began at about 0.6 s after rocket launch indicating some corona-like 
discharges before the main lightning event. 

The acoustically determined (per section 4) thunder sources are displayed 
in figure 5.4. This indicates that the lightning discharge turned south and 
extended about 4 km at a height of about 2 km. This can be compared to the 
whole-sky video-tape photograph in figure 5.5B obtained at the balloon hanger 
(fig. 1.1) about 700 m south of the rocket launcher. This indicates the main 
streamer rising and propagating to the south near cloud-base level at a height 
of about 1 km. Figure 5.6B shows the same event at about 160 ms later. Note 
in the whole-sky photographs the TV camera in the center since it is above a 
paraboloidal mirror. 

The d,(p contours from the nine pulses identified for the derivative wave- 
form in figure 5.1 are displayed in figure 5.5A. This shows a remarkably tight 
set of intersections in which all 9 contours pass through nearly the same point. 
The estimated angles from the Kiva to the source are (9,cJ)) - (30°, 177°). This 
is in quite good agreement with the estimated rocket position. The major dif- 
ference is in azimuth (177° vs 155°), but there are uncertainties in the exact 
position of the rocket because it impacted to the south near the comet observa- 
tory (fig. 1.1), despite its near vertical aiming. The slant range is estimated 
as r - 700 m based on 0 = 30° and a horizontal distance of 350 m. The same con- 
tours, except now based on the field (numerically integrated) waveform are 
shown in figure 5.6A. Here the intersection is not so tightly clustered indi- 
cating that more accuracy is generally obtained from the derivative-format data. 
This is not surprising since the analysis is based on a high-frequency approxi- 
mation and because the integration has some error. For comparison the estimated 
source location (based on the derivative data) is superimposed on the plot of 
acoustic source locations in figure 5.4. 

Now we reconstruct the source. In figure 5.7 the derivative waveform with 
the estimated r (as discussed) is used to find the peak values of the 3T/3t 
pulses. Figure 5.7B shows the individual vector pulses (as seen from the Kiva) 
oriented to the right and slightly downward. Assuming that the streamer veloc- 
ity v e ff is parallel to J and in the same direction for this case we have the 
effective reconstruction of the positive streamer in figure 5.7A. For this 
type of assumed positive streamer, denoted by v e ff 11 J we take the individual 
vector pulses in their order of occurrence and lay them out from the graph 
origin sequentially in a cumulative manner (tail of vector coincident with head 
of previous vector) . Figure 5.8 does the same thing for the T pulses from the 
field waveforms and gives a comparable picture. 

While this lightning event was a positive-going streamer, it is interesting 
to note that the impulsive nature of the fields is similar to that of natural 
negative step leaders. 
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B. Midrange Leader (s) 


Our second example is given in figures 5.9 through 5.21. As will be seen 
this is a case of significant branching. For sake of a label let us refer to 
this case as "midrange leader (s)." Figures 5.9 and 5.10 show the derivative 
fields and fields for the 20 ys record indicating leader-like behavior. 

Figure 5.11 shows the slow electric field and thunder microphone records 
from which a horizontal range of 400 m is estimated. Figure 5.12 gives the 
acoustic source reconstruction indicating an interesting streamer-path structure. 
Extending from about 3 km height, about 3 km to the east (slightly north of 
east), the pattern divides at about 2.5 km above the Kiva and descends to the 
mountain in at least three major channels. On this is superimposed three EM 
source locations determined from the derivative waveform. 

The whole-sky video-tape photographs in figures 5.13B and 5.14B show these 
three channels quite clearly as seen from the Kiva. Figure 5.14B corresponds 
to 32 ms later than figure 5.13b. 

Looking at the 0,<J> contours in figure 5.13A from the derivative waveform 
in figure 5.9, we do not see all the contours passing near a common (0,cj)). It 
appears that several sources are located above the Kiva. From figure 5.13A 
three (9,(f>) values were found as indicated on the figure. These angles are 
combined with the horizontal range of 400 m to give the estimated EM source 
locations in figure 5.12. As an aid in identifying the three sets of pulses 
corresponding to the^three locations, the directions of B 7 ? / 3 1 were considered so 
as to have all the 9T/9t pulses in a particular set have the same general direc- 
tion (as in figures 5.15 through 5.17). 

Considering the field waveform in figure 5.10 the plot of 0,<f) contours in 
figure 5.14A was made. Here four sets of pulses were identified. Sets 1, 2, 
and 3 correspond reasonably to those from the derivative waveforms as do the 
directions of the 'f pulses as compared to the corresponding 9 7 ?/9t pulses. Set 4 
in figure 5.14A was defined based on a disagreement in the direction of ^ pulses 
from those in set 1 as based on the directions in figures 5.18 and 5.21. How- 
ever, noting the small 0 values for this case the change in (j) still corresponds 
to only a small distance change on the unit sphere, but a large rotation of the 
unit vectors 1 2 » 1 3 - Identifying these set 4 T pulses with the corresponding 
9T/9t pulses it appears that set 4 pulses likely belong combined with set 1 . 

Since set 1 corresponds to small 0 an estimate of slant range based on 
400 m horizontal range is likely in error because the lightning arc is not 
exactly vertical. Based on figure 5.12 the lower heights as in figure 5.15 for 
0 = 20° are more reasonable. This is consistent with relying more on 9T/9t 
data for 0 , <J> determination. 

For this negative step leader we have v e ff J (antiparallel) . So in the 

effective streamer reconstruction in figures 5.15 through 5.21 the vectors for 
the sequential pulses are started from the graph origin with each vector head 
placed at the previous vector tail. 
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Note that values of 'f pulses depend on choice of r. For sets 1 and 4 r is 
clearly overestimated indicating the true T values are somewhat less (by per- 
haps a factor of 2) than indicated in figures 5.18 and 5.21. 

C. Nearby Leader 

Our third example is given in figures 5.22 through 5.31. This is labelled 
"nearby leader." Figures 5.22 and 5.23 show the derivative fields and fields 
for the 20 ys record. 

Figure 5.24 shows the slow electric field and thunder microphone records, 
from which a horizontal range of 95 m is estimated. Figure 5.25 gives the 
acoustic source reconstruction indicating one principal streamer channel near 
the Kiva. Extending from about 5 km height there is some branching also col- 
lecting charge about 1.5 km a little north of west. On this is superimposed 
the location of the two EM sources which were very close together on this scale. 
Figures 5.26B and 5.27B show the same whole-sky video-tape photograph indicating 
a channel to the west. Noting that this whole-sky camera is not exactly at the 
Kiva, but about 50 m north of it on the highest point of South Baldy peak (see 
fig. 1.1) the several channel location techniques give quite good agreement. 

Looking at the 0,cj) contours for derivative data in figure 5.26A we find 
two approximate sets of contour intersections. This was aided by noting the 
direction of 9T/8t pulses in figures 5.28 and 5.29. Figure 5.27A does the same 
for the field waveform pulses with good agreement, including the directions of 
the t 1 pulses in figures 5.30 and 5.31. 

Noting the two sources approximately above and below each other, these 
appear to be on or near the same main channel. It is conceivable that the 
lower source is a positive streamer and the upper streamer a negative one with 
the two streamers heading toward each other. For the plots of the_^9T/9t and T 
pulses in figures 5.28 through 5.31 we have chosen the convention v e ff J 
since such an above possibility is difficult to confirm. 


D. Return Stroke 

Our fourth example is given in figures 5.32 through 5.39. This is labelled 
"return stroke" based on the wide (1 ys) pulses as indicated in figure 5.33, as 
well as the sustained deviation of the later-time portions of the waveform. 
Figures 5.32 and 5.33 show the derivative fields and fields for the 20 ys 
records . 

Figure 5.34 shows the slow electric field and thunder microphone records, 
from which a horizontal range of 1775 m is estimated. Figure 5.35 gives the 
acoustic source reconstruction indicating two possible closures to ground. 
Extending from about 4 km height a little north of east of the Kiva the channel 
approaches the Kiva striking the ground about 2 km due east. On this is super- 
imposed the estimated locations of the EM source, one each for the derivative 
waveform and the field waveform. Figures 5.36 and 5.37 show the same whole-sky 
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video-tape photograph indicating a channel to the east, in agreement with the 
acoustic source location. 

Looking at the 0,(f) contours for derivative data in figure 5.36A we find a 
single approximate contour intersection somewhat north of east and at a high 
elevation. Similar results are found from the field waveform pulses in figure 
5.37A. These two (0,<J>) values are combined with the horizontal range to give 
the two estimates for the source location in figure 3.35. Note, however, that 
the lightning channel is not vertical and the true source location is likely to 
be farther to the east, and at a high elevation, unless some branching occurs 
back toward the west. 

^ Figures 5.38 and 5.39 give the rough streamer reconstruction based on ^ 
3T/8t and 'f respectively. The streamer is assumed positive, i.e., v e ff if J, 
in keeping with the character of a return stroke. Note, however, the high ele- 
vation. This return-stroke corresponds to a large and wide current pulse in 
the clouds. A larger range (as is possible) would give larger 3 7 f/ 3 1 and 'f' 
pulses than indicated in figures 5.38 and 5.39. 


E. Some Comments 

These four examples give representative characteristics of lightning near 
South Baldy peak. Many more recordings were obtained. It is anticipated that 
a significant fraction of these will be analyzed for future publication. These 
include both leader-like examples and return-stroke-like examples, including 
events near the horizon which appear to be more classical return strokes. 

->■ 

From the examples given the values of T have magnitudes (from transverse 
components) of about 2 x 10 12 Am/ s for return strokes to about 5 x 10 11 Am/s 
for leader-like sources (noting order of magnitudes and directions of streamers 
toward (or away from) the observer at the Kiva). Assuming a v e ff of about 
10 8 m/s gives 20 kA for return strokes and 5 kA for leader-like pulses. These 
are very variable and the number is quite approximate. Further analysis of the 
data may shed more insight into such results. 

For these data all the lightning events represent negative charge being, in 
effect, lowered to earth. It would be interesting to see what would occur in 
the rare event of positive charge being lowered to earth. 

In the recordings for the three components of field time derivative at the 
Kiva the peaks did not occur at exactly the same time. The separation of the 
sensors on the ground plane was about 11ns, and the difference in cable transit 
times was about 14 ns, giving a maximum time shift of about 25 ns. This 
resulted in peaks on the digital printout having a relative time displacement of 
up to three time intervals, or 30 ns, as was observed. For the present analysis 
the actual peaks of the individual components, regardless of time, were used. 


Unfortunately the trigger signal for the Biomation recorders was not 
recorded on magnetic tape in 1979, whereas it had been so recorded in 1978. 
This would have given even more definitiveness to the interpretations of the 
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3-component 1979 data. However, analysis of 1978 waveforms does indicate that 
to the 1 ms tape resolution, the waveforms all correspond to the first return 
stroke or the leader-like anticipation of this return stroke. 


VI . SUMMARY 


As the foregoing data and associated data analysis has indicated, fast 
nuclear-EMP related sensors and appropriate fast transient recorders have much 
to offer for the measurement and associated understanding of lightning electro- 
magnetic environments. The lightning and nuclear EMP phenomena have some simi- 
lar characteristics because of their common nature as transient electromagnetic 
phenomena. 

In this report we have discussed some of the results obtained concerning 
lightning environments from our measurements using such instrumentation on 
South Baldy peak. Three electromagnetic field components (one electric, two 
magnetic) were measured on a ground plane with 10 ns resolution for recording 
times of about 20 ys. By comparing the three field components to each other, 
and by the use of acoustic ranging, the lightning electromagnetic sources were 
approximately located and their characteristics were studied. 

The individual lightning waveforms yielded information concerning the tem- 
poral characteristics of leader electromagnetic fields. Pulse widths of the 
electromagnetic field individual pulses varied from about half a microsecond to 
less than 100 ns. Characteristic times for the rise (peak field divided by 
peak derivative) varied from something approaching 100 ns to something less 
than 30 ns; rise times were a little longer than these characteristic times for 
the rise. 

Information was also obtained concerning the temporal characteristics of 
return-stroke electromagnetic fields. Pulse widths (of the magnetic field) of 
the order of 1 ys were observed, and rise times less than 100 ns were also 
observed. Characteristic times for the rise of the order of 50 ns were 
observed . 

Comparing the fast electromagnetic-field pulses from leaders and return 
strokes at comparable distances from source to observer shows that the return- 
stroke field waveforms have generally larger amplitudes than those for leaders. 
However, for the time derivatives of the electromagnetic fields, the amplitudes 
for leaders and return strokes are quite comparable. 

Combining the three electromagnetic-field time-derivative waveforms (or 
their time integrals) one can obtain for each pulse a relation between 0 and $ 
giving a contour in the 9,<J) "plane." Intersections of such contours can give 
approximate (9,cj)) values for source location providing the corresponding pulses 
originate from nearly the same locations. With location, source orientation 
(current density transverse polarization) can then be determined. Acoustic data 
can be used to estimate distances to these sources. The source vectors l t • T 
(their peaks or peak time derivatives) can then be plotted to see how the current 


58 


direction changes from pulse to pulse. Connecting these vectors end to end 
gives a crude picture of the tortuousity of the leader path. From the magni- 
tudes of the components of T one can obtain some information concerning cur- 
rents and streamer speeds by giving these components the form Iv e f f . Typical 
current values of less than to greater than 10 kA are consistent with streamer 
speeds of the order of 10 8 m/s for leader strokes and return strokes respectively. 

This data has some implications concerning the formulation of criterion 
lightning electromagnetic environment (s) . The pulses encountered are signifi- 
cantly faster than a microsecond in their rise characteristics, even less than 
100 ns. In terms of frequency spectrum this means much more high frequency 
content for the individual pulses. A detailed frequency spectral analysis of 
these kinds of pulses would be quite useful for the formulation of lightning 
electromagnetic criteria. One should note that these waveforms are not strictly 
applicable for direct-strike lightning since measurements were not made in such 
an environment. One may expect some of the presently observed environmental 
characteristics to apply in the direct-strike case, but this is not a full 
description. 

The measurements reported here were made on South Baldy peak. Measurements 
at other geographical locations with different atmospheric conditions might 
reveal some differences in the electromagnetic waveforms from those obtained 
here. In the context of South Baldy peak more analysis of the existing data 
may reveal additional interesting features; additional waveforms are planned to 
be included in a future report. Future measurements may also obtain more com- 
plete lightning waveform information by finer time resolution, longer recording 
time windows, and/or additional and/or more accurate tie-ins with other physical 
phenomena . 
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Iron Kiva altitude: 3275m. m.s.l. 

map contour decrements: 13m 


Figure 1.1.- Top view of electromagnetic and related measurement layout 

around South Baldy Peak, New Mexico. 
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Figure 2.1.- Leader example: 20 ys scale. 
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Figure 2.2.- Same leader example: 7 ys scale. 
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Figure 2.3.- Return-stroke example: 20 ys scale 



Figure 2.4.- Same return-stroke example: 7 ys scale. 



A. Coordinate systems 
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Figure 3.1.- Coordinates for lightning data analysis. 
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Figure 3.2.- Polar 


plot for determining direction to source. 
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Date: 79228 M.S.T.: 1545:12 


Figure 5.1.- Derivative fields from rocket- triggered lightning. 
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Figure 5.2.- Fields from rocket-triggered lightning. 



Horizontal range from Kiva to rocket launcher : 350 m 


Date: 79228 M.S.T.: 1345:12 

Figure 5.3.- Slow electric field change 
and RF power received at 34 MHz from 
rocket-triggered lightning. 
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Figure 5.4.- Acoustic location of rocket- triggered lightning. 
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Figure 5.5.- 0,<t> contours for rocket- triggered derivative waveform 
and whole-sky videotape photograph. 
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Figure 5.6.- 
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B. Whole-sky photograph from hanpar 

9,()> contours for rocket-triggered waveform and 
whole-sky videotape photograph. 
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5.7.- 9T/3t for rocket-triggered lightning. 
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A. Effective reconstruction of positive streamer 
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Figure 5.8.- T for rocket- triggered lightning. 





Figure 5.9.- Derivative fields from midrange leader (s) . 
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Figure 5.10.- Fields from midrange leader (s). 
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Date: 79219 M.S.T.: 1022:51 

Figure 5.11.- Slow E field change and thunder microphone record 

of midrange leader. 
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ticks on axes at 1 km intervals. 

Date : 79219 M.S.T. : 1022:51 


Note: Three lightning channels to earth 

1,2,3 indicate set numbers of 3?/3t pulses. 


Figure 5.12.- Acoustic location of midrange leader(s). 



Figure 5.13.- 9 ,<J) 


contours for midrange leader derivative waveform and 
whole-sky videotape photograph. 




Figure 5.14.- 0 contours for midrange leader waveform and 

whole-sky videotape photograph. 
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Figure 5.15.- 9T/3t for midrange leader (set 1). 
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Figure 5.16.— 
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A. Effective reconstruction of negative streamer 
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Figure 5.17.- 3T/3t for 
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A. Effective reconstruction of negative streamer 
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Figure 5.18.- T 


for midrange leader (set 1) 
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Figure 5.19.- T for midrange leader (set 2) 
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Figure 5.20.- T for midrange leader (set 3) 
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Figure 5.21.- T for midrange leader (set 4). 
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Figure 5.22.- Derivative fields from nearby leader. 
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Figure 5.23.- 


Fields from nearby leader. 
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trigger 


Date: 79225 M.S.T.: 1246:19 

Figure 5.24.- Slow electric field change and thunder microphone record 

from nearby leader streamer. 
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ticks on axes at 1 km intervals 
Date : 79225 M.S.T. : 124619 

Figure 5.25.- Acoustic location of lightning with nearby leader streamer. 
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Figure 5.26.- 9 ,<() contours for nearby leader derivative waveform 
and whole-sky videotape photograph. 
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Figure 5.27.- 9 ,<j> contours for nearby leader waveform and 
whole-sky videotape photograph. 
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Figure 5.28. 
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A. Effective reconstruction of negative streamer 
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- 3T/3t for nearby leader (set 1) 
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Figure 5.29.- 3T/3t for nearby leader (set 2) 
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Figure 5.30.- T 


for nearby leader (set 1) 
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Figure 5.31 


T for nearby leader (set 2) 
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Figure 5.32,- Derivative fields from distant return-stroke. 



Figure 5.33.- Fields from distant return-stroke. 
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Figure 5.34.- Slow E field change and thunder microphone record 

of distant return-stroke. 



ticks on axes at 1 km intervals 
Date : 79219 M.S.T. : 1013:04 

Figure 5.35.- Acoustic location of lightning 
with a distant return stroke. 
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Figure 5.36.- 0,<t> contours for distant return-stroke derivative waveform 

and whole-sky videotape photograph. 




Figure 5.37.- 0 , cj> contours for distant return-stroke waveform and 
whole-sky videotape photograph. 
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A. Effective reconstruction of positive streamer 



B. Peaks of af/3t 
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Figure 5.38.- 3T/3t for distant return stroke. 
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Figure 5.39.- T for distant return stroke. 
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